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Abs  trac  t 


This  paper  addresses  the  problem  of  a  single  military  transport 


aircraft  that  must  carry  cargo  from  an  origin  base  to  a 


destination  base.  Refueling  is  required  and  performed  by  a 


/'tanker^  aircraft  that  originates  from  and  returns  to  a  third 


base.  The  objective  is  to  determine  the  initial  fuel  required  by 


each  aircraft  and  the  location  of  the  refueling  point  so  as  to 


minimize  the  total  fuel  consumed  subject  to  restrictions  on  the 


range  of  the  transport  and  the  tanker.  Based  on  U.S.  Air  Force 


data,  analytical  relationships  are  derived  which  allow  the 


problem  to  be  formulated  similar  to  a  constratined  spherical 


Weber  problem  with  two  main  differences:  1  )  the  objective 


function  is  non-linear,  and  <B)  some  of  the  constraints  are  a 


function  of  the  decision  variables.  Spherical  convexity  of  both 


the  objective  function  and  feasible  region  is  shown  and  used  to 


develop  an  optimal  algorithm.  Computational  experience  is  given. 
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The  problem  addressed  in  this  paper  is  that  of  a  military 


transport  aircraft  (A/C)  required  to  transport  a  cargo  of  w  pounds 
from  an  origin  base  "0"  to  a  destination  base  ”D".  The  distance 
between  these  bases  (Dao)  is  greater  than  the  range  of  the  aircraft 
thus  requiring  it  to  refuel  enroute.  Refueling  is  performed  in  mid¬ 
air  using  aircraft  called  "tankers".  The  operation  consists  of  the 
tanker  flying  from  its  base  "B"  to  some  point  along  the  route  of  the 
transport  A/C  where  it  carries  out  the  aerial  refueling.  Once  I 

refueling  is  completed  the  tanker  returns  to  "B"  while  the  transport 
continues  on  to  its  destination.  The  objective  is  to  determine  the 
initial  fuel  required  by  each  A/C  and  the  location  of  the  refueling 
point  so  as  to  minimize  the  total  fuel  consumed  by  both  A/C  subject 
to  restrictions  on  the  range  of  the  transport  and  the  tanker.  Since 
the  distances  over  which  the  aircraft  would  normally  fly  in  a  real 
application  are  large,  it  is  appropriate  that  the  geometry  of  the 
solution  approach  be  that  of  a  location  problem  on  the  sphere 
( ear  th  )  . 

The  closest  related  problem  found  in  the  literature  is  the 
spherical  Weber  problem  with  maximum  distance  constraints  formulated 
by  Ayken  (1983).  The  problem  consists  of  determining  the  optimal 
location  of  a  source  so  as  to  minimize  the  sum  of  the  weighted 
distance  from  the  source  to  a  finite  number  of  demand  points  while 
keeping  the  source  witfnin  a  distance  S.  of  a  given  point  U .  U,  may 
be  any  of  the  demand  points  and  all  points  are  locatect.on  the  surface 

of  a  sphere  (earth).  '  • 

.  *.  *-  1 
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The  optimal  solution  procedure  developed  by  Ayken  was  based  on 
properties  obtained  by  himself  and  others  for  the  unconstrained 
version  of  the  spherical  Weber  problem  (Aly,  Kay,  and  Litwhiler  1979, 


Drezner  1981.  Drezner  and  Wesolowsky  1978.  Katz  and  Cooper  1980, 
Litwhiler  and  Aly  1979).  The  problem  Me  are  considering  differs  from 
the  one  solved  by  Ayken  in  two  major  considerations.  First,  the 
objective  function  of  the  problem  addressed  in  this  paper  is  shown  to 
be  a  non-linear  function  of  the  spherical  distance  whereas  Ayken’s 
problem  dealt  with  a  linear  objective  function.  Second,  some  of  the 
constraints  of  the?  problem  we  are  considering  depend  on  some  of  the 
decision  variables.  These  properties  as  well  as  other  mathematical 
lationships  between  fuel  consumption,  distance  traversed,  cargo 
weight  and  initial  fuel  of  the  transport  A/C  are  derived  in  the 
Append i x . 

In  a  re?cent  series  of  papers,  Mehrez,  Stern,  and  Ronen  (1983), 
Mehrez  and  Stern  (1985),  and  Melkman,  Stern,  and  Mehrez  (1986),  a 
related  problem  is  studied.  They  study  the  problem  of  refueling 
tankers  escorting  smaller  aircraft  with  the  objective  of  maximizing 
the  the  effective  range  of  the  smaller  aircraft.  A  central 
assumption  in  these  papers  is  that  fuel  flow  is  a  constant.  For 
larger  cargo  aircraft,  this  assumption  is  difficult  to  support. 

Finally  it  should  be  mentioned  that  the  initial  motivation  for 
studying  the  problem  addressed,  herein  resulted  from  various 
conversations  between  the  authors  and  personnel  at  the  U.  S.  Air 
Force  Logistics  Management  Center  at  Gunter,  AFS,  Alabama. 
Subsequently,  the  A,i  r  Force  provided  raw  data,  describing  A/C 
performance  at  different  gross  weights,  a  1 t i tudes'  and  speeds  which 
was  then  transformed  into  the  mathematical  relationships'  previously 
mentioned.  Information  on  these  and  other  ldetai Is  related  to  the 
practical  origin  of  this  problem  is  also  included  in  the  Appendix. 
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3.  The  earth  is  a  perfect  sphere. 

4.  The  A/C  fallows  the  great  circle  route  (shortest  path)  in  flying 
from  one  point  to  another. 

5.  Weather  at  altitude  (i.e.,  the  jet  stream)  is  assumed  to  be 
neg 1 i g i b 1 e . 

6.  Aerial  refueling  takes  a  negligible  amount  of  time.  Thus,  the 
region  in  which  fuel  transfer  takes  place  is  considered  to  be  a 
single  point. 

7.  Everything  takes  place  at  altitude.  Thus,  take-off  and  its  cost 
are  l gnor ed . 


Notation: 

EW  =■  Transport  A/C  Empty  Weight 

MTQW  =  Maximum  Take-Off  Weight  of  the  transport  A/C 

-  Maximum  fuel  capacity  of  the  transport  A/C 

=  Maximum  fuel  capacity  of  the  tanker  A/C 

c  I  ° 

g  =  Initial  fuel  of  the  transport  A/C  (a  decision  variable) 

h  -  Initial  fuel  of  the  tanker  A/C  (a  decision  variable) 


Cargo  weight 


GW 


GW, 


Q  IN  #>l  (w) 


G 


in 


(  W  )  = 


MPF(GW) 


dr(Q ,d) 


de<B,d> 


d^( Q . d  ) 


R(g,w) 


FC ( g , w , d ) 


FN  <  w  ,  d  ) 


R“ <  h , w ) 


Gross  weight  of  the  transport  A/C  t GW  =  EW+g+w) 

Maximum  allowed  gross  weight  of  the  transport  A/C 
The  maximum  fuel  the  transport  A/C  can  carry  in  order 
to  be  able  to  take  off  with  a  cargo  of  weight  w 
MinCF  MTOW-EW-w  > 

The  maximum  fuel  the  transport  A/C  can  carry  when  it  is 
in  the  air,  given  that  its  cargo  weight  is  w 
M 1  n(  F ,,  .  GW„,„  ,,-EW-w  } 

=  Distance  traveled  by  an  A/C  in  miles  per  1,000  lbs  of 
fuel  at  a  given  gross  weight  GW 
=  The  spherical  distance  between  the  origin  base 
(8 ...  ,d,.,  >  and  the  refueling  point  (8,<z5) 

=  The  spherical  distance  between  the  refueling  point 
<0,d>  and  the  destination  (  B,., ,  do  ) 

=  The  spherical  distance  between  the  tanker  base 
<8r.  ,dr.)  and  the  refueling  point 
=  Range  of  the  transport  A/C  when  its  initial  fuel  is  g 
and  its  cargo  weight  is  w 

=  Fuel  consumed  by  the  transport  A/C  when  it  flies  a 
distance  d  and  its  initial  fuel  is  g  and  cargo  weight 
is  w 

=  Exact  amount  of  fuel  needed  by  the  transport  A/C  to 
fly  a  distance  d  when  the  cargo  weighs  W 
Cd  must  be  <_  .0 )  ] 

=  Range  of  the  tanker  when  its  initial  fuel  is  h  and 

* 

cargo  weight  is  w  (usually  w  =  0) 
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FC"'(h,w.d)=  Fuel  consumed  by  the  tanker  when  it  flies  a  distance 


d»  its  initial  fuel  is  h,  and  cargo  weight  is  w 
Cd  must  be  jC  R**(h,w>,  and  usually  w  =  03 
FN~(w,d)  =  Exact  amount  of  fuel  needed  by  the  tanker  to  fly  a 
distance  d  when  its  cargo  weight  is  w 


Cd 

must  be  _<  R~<H 

and  usually 

w  = 

0  3 

fiO.aO  -  The 

fuel  consumed 

by  the  transport  A/C  in  flying 

from 

( e,. 

» ri„,. )  to  <  0  ,  d  ) 

Cf i(0P6)  =  F  C  (  g  *  w  ,  d 

■  )  3 

f  r~  (  0  •  d  )  =  The 

fuel  needed  by  the  transport  A/C 

to 

go  f r am  ( 

0 .  d  ) 

to 

(80,60)  Cf,-<0. 

d  )  =  FN (w,dH) 3 

f,a  (  0  <  d  )  =  The 

fuel  consumed 

by  the  tanker  in  going 

from  ( 0„. 

1  ,  d  e,  ) 

to 

(  0  ,  d  )  C  f  -a  (  0  ,  d  ) 

=  F  C  "  (h,0,d;i)  3 

f ^ ( 0 . d )  =  The 

fuel  needed  by  the  tanker  to  go 

from 

(  9 ,  d  )  back 

to 

(  0i » « )  C  (  0  « 

d  )  =  FN  -  (  0  .  d  >  3 

PROBLEM  FORMULATION 

The  problem  can 

now  be  stated 

mathematical ly  as 

f 0 1 1 d ws : 

Find  <  9 . d ) 

»  g  and  h  that 

will 

(P)  Minimize 

v(0,d >  =  y~  f 
1^1 

,  (9,d  ) 

S.T.  d,  <  0 ,  d  ) 

<_  R  C  g,„... .,  (  w  )  ,  w  3 

(  1  ) 

d ,.<0.d  ) 

_  R[G,„«,,(w)  .  w  3 

(  2  ) 

d  .(  9 . d  ) 

£  R“  {  ..  ,  0  > 

(3) 

FN(w,di  ) 

<_  g  <_  g,.,,. ..  (  w  ) 

(  4  ) 

FN*  iO,0d-., 

>  y_  h  H,.,„ 

* 

(5) 

0 

<B<_  2TT 

(7) 

0 

<  6  <  TT 

<  8) 

The  first  three  constraints  represent  conditions  that  must  be 
satisfied  by  the  refueling  point.  Inequality  (1)  says  that  the 

distance  between  the  origin  base  and  the  refueling  point  should  not 
exceed  the  maximum  range  of  the  transport  A/C.  Inequality  (2)  says 
that  the  transport  A/C  must  be  able  to  reach  its  destination  from  the 
refueling  point.  Inequality  (3)  says  that  the  tanker  should  be  able 
to  make  a  round  trip  safely.  Constraints  (4),  (5)  and  (6)  are  the 
conditions  that  must  be  satisfied  by  g  and  h>  that  is,  g  must  be  at 
least  equal  to  the  amount  of  fuel  needed  to  fly  from  the  origin  to 
the  refueling  point.  Similarly,  h  should  permit  the  tanker  to  make  a 
trip  from  its  base  to  (0,d)  and  back.  Also,  the  total  initial  fuel 
for  both  A/C  must  be  at  least  equal  to  the  total  fuel  consumed. 

A  few  more  remarks  are  in  order.  First,  the  intersection  of  the 
first  three  inequality  constraints  (1),  (2)  and  (3)  represents  a 

region  within  which  all  feasible  refueling  points  must  lie.  Figure  1 
illustrates  such  a  region.  Note  that  even  though  all  feasible  points 
are  included  in  this  region,  not  all  the  points  in  this  region  are 
necessarily  feasible.  It  can  be  shown  that  any  feasible  point  must 
satisfy  (1),  (2),  (3)  and  the  following  constraint: 

g  w)-FCCg _ <w)  ,w,d  ,  3  +  CH . -FC  "  (  H,,,,. ,  Q  ,  d  •• )  ] 

£  f  ,.<  B.ri  )  +f,.  (0-,  d  )  <  9  ) 

c  -I 

wh l c h  is  derived  from  <4),  (5),  and  <6>. 
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Second,  since  values  of  g  and  h  greater  than  needed  will 


increase  fuel  consumption  unnecessarily,  we  can  restrict  attention  to 
the  case  where  the  equality  holds  for  constraint 

(6).  This  means  that  there  is  only  one  value  of  h(0,ri>  for  every 
choice  of  g (9,d>.  Consequently,  the  best  value  of  g(0,d>,  (denoted 
by  g“(9,d)>  that  will  result  in  the  lowest  total  fuel  consumption 
for  that  point  is  the  solution  to  the  following  non-linear  subproblem 
in  one  dimension  (g): 


(SP1)  Minimize  g  +  FN"  [  hr(  g  )  ,  d  3  +  h,.(g> 

S.T.  g.„lr1(0,d)  <_  g(0,ri)  <_  gm«.«<w) 


(  10) 
(11) 


where 

g,„  ri  (  9 ,  d  ) 


FN(w,d.,  )  i  f  g  r  (  B ,  eS  >  <_  0  , 

FN  (  g  r  +w  ,  d.i  )  +  g.i-  if  g  ,■  (0,s4>  _>  0  , 


(10) 


and 

g,  (8,6)  =  f«(0,*i)  +  f„<9,d>  -  CHm«„  -  FC  "  <  ,, ,  0  ,  d  ;l )  ]  (13) 

h  ( q  )  =  f  (  0  ,  b  )  +  f ,,  <  0  ,  b  )  -  C  g  ( 0 ,  d  )  -  f  ,  (  0  ,  b  )  3  (  1  7  > 

h  (  9 ,  d  )  =  F  N  "  [  h  r  (  g  )  ,  d  ]  +'  h  ,•  <  gi )  (15) 

An  important  result  due  to  Drezner  and  Wesolowsky  (1978)  says 

that  the  distance  from  a  given  point  r  within  a  circle  of  radius  TT/S 
and  center  r ,  is  a  spherically  convex  function.  Since  the  right  hand 
side  of  constraints  (1),  (2)  and  (3)  are  usually  rnugh*  less  than  TT/S, 

each  of  these  constraints  represent  a  spherically  convex  _  set.  The 

region  resulting  from  their  intersection  will  a  1  s(o  bsv:p  spherically 
convex  set.  Therefore,  Problem  (P),  without  constraint  79),  will 


have  a  spherically  convex  feasible  region. 


Our  fourth  remark  is  related  to  Aly  et.al.  (1979)  who  showed 
that  the  search  for  an  optimal  solution  to  the  spherical  Weber 
problem  can  be  restricted  to  the  spherically  convex  hull  of  the 
demand  points.  This  result  can  be  extended  to  the  case  at  hand. 

Proper  tv  1 :  The  search  for  the  optimal  refueling  point  can  be 
restricted  to  the  spherically  convex  hull  of  the  three  base  points 
(the  origin,  the  destination,  and  the  tanker  base). 

Proof :  Let  V  be  the  spherically  convex  hull  of  the  three  points  and 

X  be  a  point  on  the  sphere  such  that  X  V.  Let  d  ,  ,  d:;! ,  d be  the 

spherical  distance  between  X  and  the  origin,  the  destination,  and  the 

tanker  base,  respectively.  Aly  et  al.  (1979)  proved  that  there 

exists  a  point  X ’  £.  V  that  dominates  X  in  the  spherical  distance 

with  respect  to  the  three  base  points.  That  is,  if  d  ’  ,  d  ’  ,  and  d  ’ 

are  the  distance  between  X’  and  the  three  points,  then  d  ’  £  d,  for  i 

=  1,2,3.  Since  fuel  consumption  increases  with  distance,  Ci.e., 

FN  (  .  ,  d , , )  ,  FC  (  .  ,  .  ,  d  ,  >  «  F  N  "  (  .  ,  d  • ,  and  FC  "(.,.,  d -, )  are  all  increasing 

functions  of  the  distance],  then  the  objective  function  value  at  X’ 
is  less  than  or  equal  to  the  objective  function  value  at  X.  O.E.D. 

Taking  the  first  and  fourth  remarks  together,  it  run  be  seen 
that  the  search  for  the  optimal  refueling  point  can  be  restricted  to 
the  intersection  of  the  region  described  by  (1),  (2)'  and  (3)  with  the 

spherically  convex  hull  of  the  three  bases. 

Finally,  the  fifth  remark  concerns  the  objective  function.  It 

i  *  ** 

is  a  nonlinear  function  of  the  spherical  distance  (the  great  circle 
distance).  Using  the  fuel  ■''unction  >  see  Appendix),  it  is  seen  that 


ft  =  g  +  a/at  -  [  (a  +  a  ig)s  -  SaidJ  1  /,v  /  a 


fts  =  -a/at  +  (a'--’  +  Saide)1/S  /  a* 


f-,  =  h  +  a’/at  -  [  (a  +  a  lh)”‘  -  2a  id  ,3  1 


f  =  -a’/a,.  +  (  a  ’ +  Ea  i  d-. ) 1  /  a. 


f.  =  h+g 


-1'**  {  *c 


(  a  +  a  i  g  )a-2a  t  d  i  3  1  "'--C  a'’+2a  ,  d>-.  3  1 


(  a  +  a  ,  h  ) -  2  av  d--,  3  1  a  ’•"  +  2a  d  ,<  3 


The  way  the  objective  function  is  written,  it  is  difficult  to 
verify  its  spherical  convexity  either  by  using  the  second  derivative 
test  or  by  applying  the  definition  of  a  spherically  convex  function 


d  i  r  ec:  1 1  y 


Nevertheless,  spherical  convexity  of  the  objective 


function  can  be  shown  by  rewriting  the  objective  function  using  some 
spherically  convex  function  like  FN(.,.).  To  do  this,  define  for  any 
point  (9,0)  the  following: 


WFL  C  g  ( 9 , 6  >  ,  d  ,  3  =  weight  erf  the  fuel  left  in  the  transport  A/C  when  it 

flies  a  distance  d,  from  the  origin  to 
.•  (9,d)  given  that  its  initial  fuel  is  g(0,d). 


Thus ,  '  * 

WFLtg  (8,  it  )  ,d  ,  3  =  g(0.c5>  -  FC  [  g  (  8  .  e)  >  ,  w  .  d  ,  3 
WFL  (...)  is  spherically  convex  because  FC(  .  ,  .  )  is,  sphaer  i  c  a  1  1  y 

concave.  :: 


WFLm„„  tdi  )  =  Max  WFLC  g  <  8 .  d  )  ,  d  i  3  =  WFLC  g„,„„  (  w )  ,  d  t  3  <17) 

Note  that 

0  <.  WFLCg,.,  (0,ri  >  ,di  3  £  WFLm.„<d1),  (IB) 

and  that  WFLm*M(.)  is  spherically  convex  also. 

Since  g“<8,d)  is  the  best  value  of  g(B,d),  then  using  the 
spherically  convex  functions  FN  AND  WFL,  the  elements  of  the 
objective  function  can  be  rewritten  as 


f  !  (9,aS ) 

=  FNC  w+WFL  Cg  “'<  8  ,  d  )  ,  d  1 3  ,d  ,.  > 

(19) 

fK(0,d) 

=  FN(w»di~) 

(20) 

f «,  ( 0 ,  d  ) 

=  FN  *  <  0  ,  d ) 

(21  ) 

f  as  (  8 ,  d  ) 

=  FN  “  <  f  a+f -WFL  C  g  *  ( 0 ,  d  )  ,  d  ,  3  ,  d  » > 

(22) 

and  the  total  fuel  consumption  =  V(8,d)  =  fi(0,ei>. 

i  =  1 

Both  FN< . , . )  and  FN“ ( . , . )  are  convex  and  increasing  functions  of 
their  arguments,  which,  in  turn,  are  spherically  convex  functions. 
Therefore,  the  resulting  objective  function  is  spherically  convex 
(Greenberg  and  Pierskal la  1971,  Rockafel lar  1970). 

MOTIVATION  FOR  THE  SOLUTION  PROCEDURE  OF  ( PO > 

For  this  problem,  the  optimal  point  (8",<zS“*)  and  its  optimal 
g" <B“ ,d- )  (and  consequently,  its  optimal  h“(0-,b~))  are  to  be  found. 
Since  g“<8‘*,«5“)  depends  upon  <8,d)  and  we  do  not  know  this  in  advance 
(note  that  ( Sp 1  )  needs  to  be  solved  to  find  this  out!),  it  is 

difficult  to  use  (19),  (20),  (21)  and  (22)  directly.  ;That  is,  we  can 

not  search  for  both  (S~,e)")  and  g*(0*,cj")  simultaneously.  However, 
this  difficulty  can  be  overcome  by  creating  an  upper  bound- function 
that  can  be  improved  from  one  iteration  to  another  until-  i  t-s  minimum 
value  coincides  with  the  minimum  of  the  original  function.  To 


illustrate  this,  let 


a<8,d>  =  WFLCgO.d)  ,d,]  /  WFL„,*.,<di> 


(S3) 


then ? 

0  £  a(0,d)  £  1  for  all  <0,d> 

Thus?  g*(8,d>  would  result  in  a~(0?d)?  where  a~(0?d)  is  the  best 
value  of  a(8?d)?  i.e.? 

a~(0,d)  =  WFL C  g  ~  (  0  .  d  )  ?  d  t  3  /  WFL„.«.„(di> 

Also?  from  (25)?  we  have 

WFLCg  (0?*5 )  ?di  1  =  a  (  0 ,  d  )  WFL,„„M  (  d  i  ) 

The  elements  of  the  objective  function  can  be  rewritten  using 
a“*(0?d  )WFL,„«,„  <  d  ,  )  instead  of  WFL  C  g  "  (  0  ?  d  )  .  cl  t  3  as 
f»(0,o()  =  FNCw  +  a*<G.d)WFL,„.M(di  >  ?d  »  3 
fc(8,d>  =  FN(w?d,v) 
f,?(9.d>  -  FN  *  (  0 ,  d,v?) 

fr?(0.d)  =  FN  “  [  f  K-+f  <*-a  “  <  0  « d  )  WFL.„«.  ,•  ( d  ,  )  ?  d  -?] 

Now  setting  a(9.d)  =  a  for  all  points?  and  denoting  the  resulting 
function  by  UB(0?d?a>?  one  gets 

f  i  '  (  0  ?  of )  =  FNtw+aWFL,,,.,,.  (d  ,  )  ,d  ,  3 
f  i~’(9?d)  =  f?..<0.d>  =  FN(w,dr> 
f<?’(8,d)  =  f ...  (  0  .  d  >  s  FN“(0,d.) 
f  <  0 .  d  )  =  FN  “  C  f  e?+f  .,-aWFL,„„,,  (  d  ?  )  .  d.-.3 
where  0  £  a  £  1 

and  UB ( 0  ?  d  ?  a  >  =  2lL  ft’(0,d) 

i  =  l  ’ 

UB(0?a(?a)  has  the  following  characteristics:  *  ‘  • 

1)  It  is  spherically  convex  because  it  is  the  sum  of 

spherically  convex  functions;  ■■  '  •. 

2)  UBC8.d.a3  £  V(0,d)  ••= 

In  other  words.  since  any  g  other  then  g**(0,d)  would  result  in  a 
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(24) 

(25) 

(26) 


(27) 
(20) 
(29) 
(  30) 

(  31  ) 


(32) 


higher  cost  for  V(8»d> »  and  also  produce  an  "a"  that  is  different 
from  a~(0,d)»  then  any  "a"  used  in  UB(8»d,a)  would  make  UB(0,d.a) 
<_  V(8,d).  Therefore  UB  is  an  upper-bound  for  the  objective 
function; 

3)  UBC0 , 6  ,  a“  ( 0 ,  d  >  3  =  V(0,d)  for  any  <0,«5>.  That  is,  the  equal 

sign  of  eq .  (32)  holds  true  when  a  =  a“ (0,d ) 

Therefore 

UBC0-  ,6-  ,a~  <0“  ,d~  )  3  =  V  (  0  “  ,  d  “* )  (33) 

4)  UB(B.d,a)  =  V(0,d)  at  the  boundary  points  where  di(0,d)  = 
RCgm„...  (w)  ,  w3  because,  at  those  points,  WFL„,„ (  d  >  =  0.  Thus 
the  value  of  a  does  not  matter. 

How  can  we  use  this  UB  function  to  solve  (P)?  Given  a 
particular  value  of  a,  we  can  optimize  the  UB  function  to  find 
the  best  (0,eS>  that  corresponds  to  it.  Now,  if  at  some 

iteration,  we  have  the  value  of  a"(8“’,6“)  [without  knowing 
(0",d*)3,  and  we  optimize  the  UB  function,  then  we  will  end  up 
with  (9**,d").  This  is  true  since  if  (0"»d“)  minimizes  the  true 
objective  function,  it  will  also  minimize  UBC 0 , d , a* ( Q" , d " > 3 . 

SOLUTION  PROCEDURE  OF  (P) 

1 .  Initial ization; 

Let  k=0,  and  a-'=l.  Start  with  a  point  (0'\*rio>  that  satisfies 
(1),  (2)  and  (3)  [(See  Yamani  1986  for  the  selection  of  a  good 

starting  point). 3 

2.  Let  k=K  +  l  .  Find  the  point  that  will  minimize  UB  <  0  ,"15.  ,,a‘ ~ 1  ) '  s  .  t  . 

(1),  (5),  and  (3).  Let  the  solution  be  (0'',d'r>. 


3.  For  (0‘-  ,«J‘)  ,  solve  (SPl)  to  find  the  best  g"(8'  ,d'  ).  Use 


equation  (25)  to  find  a~(  Sk"  *  >  . 


4.  If  aM8*,d*>  *  a*'"1  ,  we  have  an  improvement  in  the  objective 
function.  Therefore!  let  a1'  =  a“  ( 0*'  , d*  )  and  go  to  step  2. 
Otherwise,  the  solution  is  (0*'»d*'),  g*  <  0*r  ,  6tr  )  »  and  h~ ( 8k , d*  ) 
is  calculated  from  eqs.  (14)  and  (15). 

Convergence  Proof 
We  need  to  show  that 


2^  f  i  <0*~*  ,d*~l  )  >  ^f^GSsf1')  (34) 

i  =  1  ”  i  =  1 

Assume*  at  iteration  k-1,  that 
a*- (g*-t  >  *  aV!-«. 

(Otherwise,  we  would  have  stopped  there.) 

By  setting  a'"""11-  =  a*  ( S*-1  ,  d*"1  )  *  we  have 

f*<8 *-1,6*-*)  =  UB(e>'“1  ,a*~r  )  (35) 


At  iteration  k,  since  (8|r,dl">  minimizes  UB  (  8  ,  d  ,  a l' , 


UB  (  B’’  -  1  ,  d'  : -1  ,  a^':  “  :l  )  >  UB  (  0*'  ,  d1'-  ,  a*  -1  > 


Bu  t  , 


UB  <  01'  ,  d,,; » a1''"1 )  >  ^If  ,  (B'*;,dK )  ■ 

~  i  =  1 

and  inequality  (34)  follows. 

Moreover , 

f .  (  S' . 1  ,d"  '  >  1  f,  (0‘-  ,d1'  )  , 


(36) 


(37) 


(33) 


i  =  1 


i  =  l 
4 


if  U8(0*"  ,d*  ,a*";~- J  )  >.  f ,  <  B'<  ,  d*’ >  .  -  (39) 

The  last  inequality  holds  true  if  a  ( 0t,:  ,  f5tr  )  a1'-'"1.  Therefore,  the 

true  objective  function  improves  from  one  iteration  to  -another. 

Since  the  problem  is  spherically  convex,  the  solution  procedure  will 

converge  to  the  global  optimum  solution  (Drezner  and  Wesolowsky 
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1973).  Q.E.D. 


Numerical  Results 

The  above  procedure  was  programmed  on  a  Vax  11/750  in  BASIC.  There 
was  no  attempt  at  program  efficiency.  Several  test  problems  were 
solved  to  represent  a  variety  of  possible  geographic  configurations 
and  complexities.  The  minimum  CPU  time  was  approximately  0.5  second 
(2  iterations)  and  the  maximum  was  almost  1.0  second  (4  iterations) 
for  convergence  within  one  degree  longitude  and  latitude.  Convergence 
within  one  minute  longitude  and  latitude  required  only  one  more 
iteration  in  each  case  with  a  maximum  CPU  time  of  1.3  seconds.  The 


results  of  these  test  problems  are  summarised  in  Table  2,  with  Table 
1  providing  a  key  for  understanding  the  entries  in  Table  2. 

Summary 

In  this  paper  we  have  shown  that  the  single  aircraft  refueling 
problem  is  a  non-linear  convex  location  problem  on  the  sphere.  The 
algorithm  that  has  been  presented  is  an  optimal  solution  procedure 
for  this  problem.  In  the  sequel  to  this  paper,  these  results  are 


used  as  the  basis  for  solving  the  more  complex  multi-aircraft 


refueling  problem. 


TABLE 

1 :  Key  for  trial  runs 

chart  of  Table  2 

Location 

Longitude 

Latitude 

Code 

New  Jersey 

7  5W 

AON 

A 

Delaware 

75W 

38N 

B 

North  Carolina 

78W 

35N 

C 

Puerto  Rico 

66W 

18N 

D 

Azores  Islands 

25W 

37N 

E 

Iceland 

20W 

65N 

F 

Germany 

10E 

SON 

C 

Turkey 

30E 

AON 

H 

Saudi-'Arabia 

A7E 

25N 

I 

Egypt 

28E 

30N 

J 

England 

OE/W 

52N 

K 

W  =  Cargo  Weight 

Code 

100,000  lbs 
'200,000  lbs 

Fuel  ACCR  =  Fuel 

1 

2 

Accuracy 

10  means  to 

within  10  lbs 

POS  ACCR  *  Position  Accuracy 

'  C  •  ♦  ° 

- 

1  means  to  within  1  degree 


ble  2:  Trial  runs  for  Test  Problem 


Appendix  A 


Derivation  and  Characterization  of  the  Fuel 
Consumption  Function  and  Other  Related  Functions 

Here*  mathematical  representation  of  the  various  relation¬ 
ships  between  fuel  consumption,  distance  traveled,  cargo  weight, 
and  initial  fuel  of  a  typical  jet  A/C  are  derived. 

The  United  States  Air  Force  provided  raw  data  describing  A/C 
performance  at  different  gross  weights,  different  altitudes.  and 
different  speeds.  These  raw  data  have  to  be  transformed  before 
they  are  useful.  A  sample  of  these  data  is  shown  in  Figure  8  for 
a  C-5A  transport  A/C  flying  at  an  altitude  of  31,000  feet. 
Different  curves  are  presented  for  different  A/C  weights,  and 
they  show  the  distance  traveled  per  1,000  pounds  of  fuel  burned 
at  a  given  speed  and  at  that  given  gross  weight  (GW).  Consider 
the  points  on  the  curve  indicating  99  percent  maximum  specific 
range.  It  is  assumed,  without  loss  of  general i 1 ty ,  that  the 

transport  A/C  will  operate  at  those  points  since  speed  is 
maximized,  while  fuel  consumption  is  only  1  percent  greater  than 
the  minimum  possible.. 

The  distance  traveled  in  miles  per  1,000  pounds  of  fuel 
burned  when  the  A/C  gross  weight  is  GW,  denoted  here  as  MPF ( GW ) , 
is  plotted  against  gross  weight  GW  in  Figure  3  using  the  99 
percent  line  in  Figure  8.  Least  square  linear  -anil  quadratic  tits 
to  the  data  were  tried,  and  both  fit  well  as  seen" in  Figure  3. 

Forthelinearfit, 

.  ’  "£■-  •  » 

MPF  (GW)  =  at.  +  a ,  GW ,  '  := 

wher  e 

a.. 


36.2839,  a 


-0.087 


Air  Nautical  Miles  per  1,000  Pounds  of  Fuel 


and  correlation  coefficient 


-0.9919 


For  the  quadratic  fit, 
MPF(GW)  =  b,_.  +  b  i  GW  +  b^-GW** 

where 


b,_,  =  43.7616,  b.,.  =  -0.0576,  bG  =  E .  94  x  10  ° 
and  correlation  coefficient  =  -0.9983. 

Only  the  linear  form  of  the  function  will  be  used  here.  The 
use  of  the  quadratic  form  is  similar  ( Yamani  1986). 

Note  that  at  any  moment  during  the  flight, 

GW  =  EW  +  w  +  f 

where  f  =  the  present  (instantaneous)  amount  of  fuel  and  EW  =  A/C 
empty  weight.  During  flight,  f  changes  due  to  fuel  burning; 

thus,  GW  changes  while  EW  and  w  stay  the  same.  So, 
dGW  =  df 


The  Range  Function  R < g , w ) 

Initial  Gross  Weight  =  EW  +  w  +  g 
Final  Gross  Weight  =  EW  +  w 


Then , 


R  (  g  ,  w  )  = 


EW+ w+g 
EW+.w 

g 

o 

g 

o 


J 

l 


NPF(GW)  dGW 


NPF(GW)  df 


MPF (EW+w+f ) 


df 


The  Fuel  Consumption  Function  FC(g,w,d) 
Let  gr  equal  the  final  amount  of  fuel  left. 
Then,  '  "£•• 


FC ( g , w , d )  =  g  -  g. 


Thus,  finding  FC  reduces  to  finding  gr. 


which  can  be  accomplished 


by  solving  for  gr  in  the  following  equation: 


d  =  \ MPF(EW+w+f )  df 

gr  J 


d  =  \Ca0  +  a1(EW+f+w)3  df 

J 


which  results  in 


g  r  =  -a /a,.  +  C<a+a;ig>‘-  -  Said]1''®  /  al 


Thus , 


FC(g,w,d)  =  g  +  a/a ,  -  C  (a  +  aig)e  -  2aid])V'"  /  al 


where 


a  =  a  o  +  a  i,  <  EW+ w  ) 


The  Fuel  Requirement  Function  FN(w,d> 


If  the  distance  d  to  be  flown  and  the  cargo  weight  w  are 


known  in  advance,  then  the  amount  of  initial  fuel  g  must  be  at 


least  as  great  as  the  amount  of  fuel  needed.  If  g  is  set  equal 


to  FN<w,d)»  then  we  can  use  the  range  function  R(g,w)  to  solve 


for  g.  That  is,  set 


R(g,w)  =  d  =  la’  +  a  ,  w  +  <  a  ,  /  E  >  g  3  g  ■ 


where 


a’  =  a,,  a,EW 


Solving  for  g  results  in  the  following  equatib’n: 


FN!w,d)  =  g  =  -w-a’/a,  +  [(a’+a,w)r"  +  Sajd]1'''”  /  a 


Table  3  summarizes  the  various  fuel  functions.'1 


'.,0  .v.w/a' 


TABLE  3:  Summary  of  the  fuel-related  functions 
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